Double-layer quantum Hall antiferromagnetism at filling fraction v = 2/(odd integer) 
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A low energy action for double-layer quantum Hall systems at filling fractions u — 2/m (m an odd 
integer) is introduced. Interlayer antiferromagnetic exchange induces a phase with canted spin order, 
and also a spin-singlet phase. Universal properties of zero and finite temperature transitions are 
obtained. We compute the critical temperature at which the canted order vanishes in a Kosterlitz- 
Thouless transition. Implications for recent light scattering experiments at v — 2 are noted. 
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There has been much recent work on double-layer 
quantum Hall systems, the majority of which has fo- 
cussed on the case where the electron tunneling rate be- 
tween the two layers is small Then, the electron 
layer index plays the role of a pseudospin, and for the case 
where the total filling factor v = l/m (to an odd inte- 
ger), very interesting new physics arises from long-range 
correlations in the pseudospin orientation. However, the 
tunneling term acts like a 'magnetic field' in pseudospin 
space, and so spontaneous long-range pseudospin order, 
and the associated finite temperature (T) phase transi- 
tion, is only possible when the tunneling is vanishingly 
small §. 

Stimulated by recent light scattering experiments |^] 
at v = 2, we present here a general low energy theory for 
double-layer systems at filling v = 2/m in the presence 
of moderate interlayer tunneling. We find a rich phase 
diagram with interesting transitions both at T = and 
T > 0. In contrast to the phases at v = l/m, which 
are driven by ordering in pseudospin space, the phases 
at v — 2 /to are associated with ordering in the physical 
electronic spin space. Consequently, our order parame- 
ters are defined even in the presence of interlayer tun- 
neling; indeed, moderate interlayer tunneling is required 
to stabilize some of our phases. We will use our results 
to interpret recent experiments ||, and argue that they 
show indirect evidence for our T > phase transition. 

It is useful to begin discussion of the physics at v = 
2/m by considering the case where the layer separation, 
d, is much larger than the magnetic length, l Q . Then 
the two layers (labeled 1,2) are approximately decoupled, 
and each separately has filling fraction v\ = v<i = l/m. 
Their ground states will be the familiar Laughlin states 
for to > 1 , or a fully filled lowest Landau level at to = 1 , 
both of which have a large energy gap to all charged 
excitations Q . These states are also fully spin polarized 
and there is significant intralayer ferromagnetic exchange 
P,|]0)- The low- lying excitations in each layer are spin 
waves which have a small excitation gap given precisely 
by the Zeeman energy gpsH (the gyromagnetic ratio g 
and the Bohr magneton [Ib will henceforth be absorbed 
by a rescaling of the magnetic field H). For small g, a 



complete description (|J^] of the low energy excitations 
of each layer can be given in terms of an action for unit 
vector fields ni,2 ( n\ 2 = 1) representing the orientation 
of the ferromagnetic orders. Spin waves are small fluc- 
tuations of fi\^2 about an ordered state, while charged 
quasiparticles are skyrmion textures of n\p- 

Now reduce the value of d and couple ni and n,2 ■ The 
simplest allowed coupling between them is an antifer- 
romagnetic exchange interaction. These considerations 
lead to the following imaginary-time (t) effective action 
(in units with h = ks — 1) 
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The intralayer ferromagnetic spin correlations |]J^],|| are 
controlled by Cp: M = 1/47ttoZ 2 is the magnetiza- 
tion density per layer, p° s is the spin stiffness of each 
layer when they are well separated (for m = 1, we 



have M p° s — e /(16v27reZ ) and A accounts for the 
Berry phase accumulated under time evolution of the 
spins (eijkdAk(n)/dnj — ni). The interlayer antiferro- 
magnetic correlations are induced by the positive cou- 
pling J ~ M A 2 as /[/ where A sas is the tunneling ma- 
trix element between the layers, and U ~ e 2 /el a is the 
Coulomb interaction energy. 

Some potentially important terms have been omit- 
ted from So and our analytic computations: the Hopf 
term which endows the skyrmions with fractional statis- 
tics, and the long-range Coulomb interaction between the 
skyrmions. We believe this is permissible because of the 
charge gap noted earlier. Further fTpf , as the layers are 
antifcrromagnetically correlated, skyrmions in one layer 
will be correlated with anti-skyrmions in the other, and 
this neutralizes the leading contribution of both terms. 
This latter argument should continue to hold even if the 
charge gap were to vanish at a quantum critical point. 
Note that no new term is necessary to induce charge 
transfer between the layers: a hedgehog/anti-hedgehog 
pair in the two layers corresponds to an event transferring 
skyrmion number between them. Such spacetime singu- 
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larities are absent in the strict continuum limit but ap- 
pear when a short-distance regularization is introduced. 
Finally, for m > 1 and larger g, the spin Laughlin quasi- 
particles become the lowest energy charged excitations, 
but these can be neglected for similar reasons. 
Now we parameterize 



rtt = (-!)*(! -£ 2 ) 1/2 ™ + £ 
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where the constraints n\ 2 = 1 are now replaced by 

n 2 = 1 and L ■ n = 0. Because the layers are antifer- 
romagnetically correlated we expect that L will not be 
too large. We insert (||) into ([!]), expand to quadratic 
order in L, and then integrate out the L degrees of free- 
dom. This yields the following effective action for the 
antifcrromagnetic order parameter n 
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where t = (J/2p° s M§) 1 / 2 and c = (2p° s J/M§) 1 / 2 . This 
is precisely the action of the 2+1 dimensional quantum 
0(3) non-linear sigma model in a field H coupling to the 
conserved global 0(3) charge. It is expected to apply 
to double-layer quantum Hall systems with v = 2/m at 
length scales larger than A -1 ~ l a . 

For the special case m = 1, we will sharpen the 
quantitative theoretical predictions of Si by a Hartree- 
Fock (HF) analysis of a realistic, microscopic double- 
layer Hamiltonian jy]. The HF theory will be used to 
compute renormalized T — energy scales which com- 
pletely specify the correlators of Si at low T: in this 
manner we determine observables of the system with no 
free parameters (for m > 1 these energy scales remain as 
phenomenological parameters). We use the Hamiltonian 
H = H + Hi with 



Ho 



h.c. 



H 



where Ci aa annihilates an electron in the lowest Landau 
level in layer i (i = 1,2) with spin a (a = ±1) in the 
z direction (we assume H = (0,0, if)) and with intra- 
Landau level index a. Interlayer tunneling induces the 
symmetric-antisymmetric energy separation A sas . The 
Coulomb interaction part of TL is 



g q 2 l 2 /2 ei q m ( ai - a2 )l 2 



fit ft f< 
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where q is a wavevector, f2 is the area of the sample, and 
the interaction potentials are Vij = 2-7re 2 /eg for i = j and 
Vij = (2ire 2 /eq)e-i d for i^j. 

The T — phase diagram (l^Jl^l for the action S\ is 
shown in Fig. |l|. For v = 2& topologically identical phase 
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FIG. 1. T — phase diagram of Si . The phases are pictori- 
ally represented by the orientation of the spins in the two lay- 
ers, with H pointing vertically upwards; C (N) has a broken 
0(2) (0(3)) spin symmetry. There is a Kosterlitz-Thouless 
transition at T = T c > in C. The inset shows the phase 
diagram obtained from a microscopic HF calculation at v — 2 
(m = 1), where the asterisks represent the two experimen- 
tal samples of ref. We argue in the text that the HF 
theory overestimates the stability of the C phase, and that 
experiments suggest that the actual SS region encloses the 
left sample in the inset. 

diagram is obtained in the HF computation and is 
shown as an inset. The quantum phase transitions be- 
tween these phases are continuous and are accompanied 
by the softening of the intersubband spin density excita- 
tions. The phases are described below: 

(I) Fully Polarized Ferromagnet (FPF): In Si this is 
present for H >^ ctA 2 . This phase has (n\ z ) — (ri2z) — 
1. It is continuously connected to the large d limit dis- 
cussed earlier. 

(II) Canted (C): We now have (ni z ) = (n 2z ) ^ 0; ancl , 
for example, (ni x ) = —(n2x) 7^ 0. This phase has a bro- 
ken spin rotational 0(2) symmetry in the x — y plane. 
For m = 1, the HF phase boundary between the FPF 
and C phases is V_/A sas = (A sas /H)[l - (#/A sas ) 2 ], 
and that between the C and SS phases is V_/A sas = 1 — 
(ff/A sas ) 2 , where V± = £ £q e-» 2 '°/ 2 [Vn(g) ± V 12 (q)}; a 
wavefunction for the C phase is obtained by the standard 
HF methods. For m > 1 a caricature of the wavefunc- 
tion is two separate Laughlin states at v% = U2 = 1/m 
but polarized in the orientations shown. 

(III) Neel (N): This is the limiting case of C with 
{ni z ) — ( n 2z) — achieved at H = 0. Now an 0(3) 
spin rotation symmetry is broken. 

(IV) Spin Singlet (SS): This corresponds to the quan- 
tum disordered phase of the 0(3) sigma model. The 
ground state is a spin-singlet and is therefore unaffected 
by H: its wavefunction is the same as that at H = 0. 
For m — 1, in the independent electron HF picture, the 
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electrons fill the layer-symmetric subband, with spin-up 
and spin-down levels equally populated. However, it is 
well-known that HF theory overestimates the energy of 
a nonmagnetic phase like SS because correlations be- 
tween opposite spin electrons, important for reducing the 
Coulomb energy, are now absent. It is likely, therefore, 
that the real SS phase is stable over a larger parameter 
region than that in our HF approximation, but the topol- 
ogy of the HF phase diagram in the inset of Fig. [I] (which 
is identical to that for Si) should be correct. To build 
in charge correlations, one can use an approach similar 
to the Heitler-London picture of the hydrogen molecule, 
and consider pairs of electrons with their charge localized 
in opposite layers, while their spins form singlet bonds. 
Indeed, such a charge-localized picture was behind our 
introduction of the actions <So,i- In such an approach, 
an alternative wavefunction for the SS phase (valid for 
m = 1 and m > 1) can be obtained in the J — > oo limit: 
pairs of electrons in opposite layers bind to form spin sin- 
glet, charge 2e bosons, which then condense into a boson 
Laughlin state at filling fraction l/2m, as demanded by 
the strength of the magnetic flux. 

It is worth noting explicitly here that the HF compu- 
tations at m = 1 allow us to assert that all the different 
phases of S± are the ground states in realistic parameter 
regimes. For m > 1, it remains an open question as to 
whether the phases of Si other than FPF are accessible, 
although we consider it a likely possibility that at least 
C will exist. 

We now turn to the physics at T > 0. Only the N 
and C ground states have a broken spin rotation sym- 
metry; the 0(3) symmetry of the former implies that the 
symmetry is restored at any T > 0, while the 0(2) sym- 
metry of the latter implies a Kosterlitz-Thouless phase 
transition at a T = T c > 0. We may characterize the 
order parameter fluctuations in both phases by a T = 
spin stiffness p s (H) such that the energy cost of rotations 
of the order parameter by a slowly varying angle 4>{r) is 
Ej, = (p s (H)/2) J d 2 r\ V0(r)| 2 . A crude estimate ||fll) 
of T c is T c ps p s (H), although this must fail as H — ► 0. 
In the latter limit it is possible to obtain an exact lead- 
ing asymptotic result [jj T c = 2irp s (Q)/\n{p s {Q)/H) for 
ln(/) s (0)/H ) 3> 1. For m — 1 we computed p s {H) in the 
HF calculation and the results are shown in Fig. ||. We 
see that the T c estimates are well in the experimentally 
accessible regimes for typical GaAs-based semiconduc- 
tor samples. We emphasize that the Kosterlitz-Thouless 
transition at T c is present even in the presence of inter- 
layer tunneling, unlike the case for the pseudospin tran- 
sition nu a t v — 

A more precise approach to the T > properties is to 
expand in the deviation from the H = T = quan- 
tum critical point between the N and SS phases at 
t = t c ~ A -1 . This quantum critical point is described 
by a renormalizable quantum field theory (with upper- 
critical spatial dimension d = 3), and so all thermody- 
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FIG. 2. Ground state spin stiffness p s (H) of a v = 2 dou- 
ble-layer system obtained from the microscopic HF calcula- 
tions. It is non-zero only in the N and C phases. The solid 
line has H — 0, the dashed line H = 0.05e 2 /eZ o , and the dot- 
ted line H = 0.08e 2 /eZ o - The layer separation is d — 1.0l o and 
we have also included corrections from finite layer thickness 
d w — O.Slo- In this figure, V- is fixed to be 0.36e 2 /eZ o for the 
given values of d and d w , and p 3 (H) is shown as a function of 
A sas for several values of H. (In typical GaAs-based samples, 
e 2 1 d is on the order of 50 — 100A', which gives p s ~ 1 — 2K.) 

namic properties are universal functions of energy scales 
characterizing 'relevant' perturbations from this critical 
point; corrections due to irrelevant operators require ad- 
ditional energy scales and will be neglected here. Two 
of the relevant energy scales are the 'bare' couplings T 
and H (there is no renormalization of the scale of H be- 
cause it couples to a conserved charge [[l2|) and a third 
(the last) measures deviation of t from t c . For t > t c 
we choose [[uj this energy scale to be A, the energy 
gap of the SS state at T = H = 0, while for t < t c 
we choose the renormalized spin stiffness p s (0) of the N 
state also at T = H = 0. As t approaches t c we have 
A ~ (t - t c y, while p s (0) - (t c - i)^- 1 )", where v 
is the correlation length exponent of the classical three- 
dimensional 0(3) ferromagnet. For m = 1, the micro- 
scopic HF calculation gives A = A sas — V-/A sas , and 

Ps (0) = ((1 - (A sas /V^) 2 )/8^)E q (^) 2 e- 92i ^ 2 I4i((z). 
Notice that these are consistent with the mean field ex- 
ponent v — 1/2 in the upper-critical dimension d = 3. 

One of our main results, which follows from the con- 
siderations above, is that the critical temperature T c at 
which the ordering of the C phase disappears obeys, for 
t > t c , 

T c = J ff* > (A/ J ff). (4) 

Here \l/>(it) is a universal function of u with no arbitrary 
scale factors, and obeys the exact relation ^>(u > 1) = 
(because [ pdp2] ] the T = boundary of the SS phase 



3 



is given precisely by the condition A = H). A similar 
scaling form holds for t < t c with T c = < (p s (0)/ H). 
We computed the functions \&>,< in an expansion in e = 
3 — d using recently developed methods |ll| and found to 
leading order 



*>(u) = [33(1-u 2 )/(10ttV)] 



1/2. 



(5) 



the structure of the subleading terms is quite compli- 
cated and is similar to that discussed elsewhere p5| . 
This result is valid for all u, except for u very close 
to 1; in that case we find, by a mapping to the dilute 
Bose gas problem, the exact asymptotic result [ p^|Jl6| 
^>>{u — ► 1) = y ln(l/y)/(41nln(l/y)) with y = 1 - u, 
which holds for ln(l/y) 3> 1. For t < t c the e expansion 
holds for ^ < (u/^/e) and we obtained 



*<(u/Ve) = [(33 + 3 M 2 )/(10^ 2 e)] 1 / 2 



(6) 



Again this result is valid for all u, but now fails for u — > 
oo (which is H -> 0). While T C (H = 0) > for all 
e < 1, we noted earlier that T C (H = 0) = for e = 1; 
the latter property will not appear at any order in the 
e expansion. Using results special to d — 2 for H — > 
discussed earlier, we have instead the exact asymptotic 
form x E'<(u — > oo) = 2iru/ lnw. 

We draw attention to a particularly simple and strik- 
ing limit of the above results. At t = t c we have T c = JCH 
where K. = ^>(0) = , 3 / <(0) is a universal number. Fur- 
ther, we do not expect any large or singular variation in 
T c if t is close to but not exactly t c . As both H and 
T c are directly measurable energies, this relationship is 
amenable to a direct experimental test. On the theoret- 
ical side, while at present there is only the leading term 
in a e expansion for the value of JC, it should be possi- 
ble to obtain a reasonably precise result using quantum 
Monte Carlo simulations of double-layer lattice spin sys- 
tems Jl^,|l^], which have been limited to H = so far. 
Universality implies that these lattice models will have 
the same value of fC as the quantum Hall system, and it 
appears to us that the simulation for H ^ should also 
be free of the fermion sign problems. 

We have also obtained results in the e expansion for the 
crossovers of the dynamic spin susceptibility at frequency 
w as universal functions of the energy ratios lu/T, H/T 
and A/T (p s (0)/T): the methods are similar to those of 
ref. |lq| , and results will be presented elsewhere. 

We turn now to a comparison with recent light scatter- 
ing experiments ||. The high density sample (the right 
sample in the inset of Fig |l|) shows 'mode-softening' con- 
sistent with a T > phase transition which we identify 
with that above our C phase. Using input parameters 
from the HF calculation in (||), we obtain the prediction 
of T c ~ 0.5K, to be compared with the experimental 
value T c ~ 0.52K: the good agreement must be consid- 
ered fortuitous until the accuracy of the s expansion is 
better understood. The same sample also shows marked 



T dependence for T > T c in the light scattering spectrum 
at lo of order or greater than T: a natural explanation for 
this could be a crossover into the 'high-T" region above 
the quantum-critical point ]l5p . The lower density sam- 
ple (the left sample in the inset of Fig [l]) shows no mode 
softening and little T dependence in the light scattering 
spectrum: we suggest that this sample is in the SS phase. 
The HF computation puts this sample in the C phase, 
but as we discussed earlier, this could be in error because 
the HF theory overestimates the stability of the C phase. 

Finally we note that we expect similar considerations 
to apply to all double-layer quantum Hall systems with 
v = 1v\ where a single layer at filling v\ forms a fully 
polarized quantum Hall state with a charge gap. 
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